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SOME  THREE-DIMENSIONAL  INCLUSION  PROBLEMS  IN  ELASTICITY1 

by 

M.  K.  Kassir2  and  G,  C.  S 1 h 3 


Abstract 

The  theory  of  potential  functions  is  applied  to  solve  a 
number  of  three-dimensional  problems  involving  sheet-like 
inclusions  embedded  in  elastic  solids.  Two  types  of  inclu¬ 
sions  are  considered;  namely,  that  of  a  rigid  elliptical  disk 
and  a  rigid  sheet  containing  an  elliptical  hole.  By  varying 
the  ellipticity  of  the  disk  and  hole,  certain  information  on 
the  general  character  of  the  stresses  around  a  plane  inclu¬ 
sion  of  arbitrary  shape  may  be  obtained.  More  precisely,  if 
reference  is  made  to  a  suitable  coordinate  system,  the  func¬ 
tional  forms  of  the  stresses  in  the  close  neighborhood  of  the 
inclusion  border  can  be  expressed  independently  of  uncertain¬ 
ties  of  both  the  inclusion  geometry  and  of  the  applied 
stresses  or  displacements.  In  general,  the  intensification 
of  the  local  stresses  can  be  described  by  three  parameters 
which  may  be  used  to  establish  criteria  for  the  failure  of 
the  solid  containing  the  inclusions. 

]This  research  was  supported  by  the  U.S.  Navy  under  Contract 
Nonr- 61 0 ( 06 )  with  the  Office  of  Naval  Research  in  Washington, 
D.C. 

^Department  of  Civil  Engineering,  The  City  College,  New  York, 
New  York. 

^Professor  of  Mechanics,  Lehigh  University,  Bethlehem,  Penn¬ 
sylvania. 


Introduction 


During  the  past  few  decades,  considerable  attention  has 
been  devoted  to  the  solution  of  two-  and  three-dimensional 
problems  of  stress  concentrations  around  inclusions  of  a 
variety  of  shapes.  Since  the  literature  on  this  subject  is 
exhaustive,  only  those  works  which  are  pertinent  to  the 
present  study  will  be  cited. 

The  problem  of  a  thin  rigid  circular  disk  embedded  in  an 
infinite  solid  and  subjected  to  a  constant  displacement  nor¬ 
mal  to  its  plane  was  solved  by  Collins  [1].  His  results  are 
equivalent  to  the  slow  steady  motion  of  a  rigid  disk  in  a 
viscous  fluid.  In  a  recent  paper,  Keer  [2]  has  considered  a 
similar  problem  in  which  the  disk  is  displaced  in  its  own 
plane.  The  case  of  an  infinite  solid  containing  a  rigid 
sheet  with  a  circular  hole  was  also  discussed  in  [2],  The 
disturbance  of  an  ellipsoidal  inclusion  in  an  otherwise  uni¬ 
form  stress  field  was  examined  by  Eshelby  [3,4].  In  the 
limit  as  one  of  the  principal  axes  of  the  ellipsoid  vanishes, 
the  solution  to  the  problem  of  a  flat  elliptical  disk  may  be 
deduced  from  the  work  in  [3,4]. 

For  the  purpose  of  assessing  the  strength  degradation  of 

solids  due  to  the  presence  of  disk-shaped  inclusions,  it  is 

important  to  have  a  knowledge  of  the  singular  behavior  of 

the  stresses  near  the  sharp  edges  of  the  inclusions.  To  this 
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end,  the  present  investigation  is  concerned  primarily  with 
the  determination  of  stress  solutions  of  the  following  bound¬ 
ary-value  problems: 

(1)  A  plane  inclusion  of  elliptical  shape  in  an  other¬ 
wise  uniform  tensile  field. 

(2)  Elliptical  disk  displaced  in  its  own  plane. 

(3)  Displacement  given  to  a  rigid  sheet  with  an  ellip¬ 
tical  hole. 

(4)  El  1 iptical ly-shaped  disk  displaced  out  of  its  own 
pi ane. 

Referring  to  a  system  of  Cartesian  coordinates  x,y,z, 
the  z-axis  will  be  directed  normal  to  the  plane  of  disconti¬ 
nuity  which  is  bounded  by  the  ellipse 

x2/a2  +  y2/b2  =ltz-0  (1) 

where  a  and  b  are  the  major  and  minor  semi-axes  of  the 

ellipse,  respectively.  The  center  of  the  ellipse  is  located 

at  the  origin  of  the  coordinate  system.  The  rectangular 

components  of  displacement  ux,  uy,  uz  and  stress  oxx,  °yy> 

— ,  T  are  assumed  to  be  continuously  differentiable  at 

all  interior  points  of  the  solid  and  take  definite  values  on 

either  side  of  the  ellipse  except  that  on  the  periphery  of 

the  ellipse  the  stresses  may  become  infinitely  large.  At 
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large  distances  from  the  origin,  all  the  stresses  and  dis¬ 
placements  tend  to  zero.  The  problem  is  to  find  a  suitable 
solution  of  the  Navier's  equation  of  linear  elasticity  for 
a  homogeneous,  isotropic  body. 

In  the  absence  of  body  forces,  the  displacement  vector 
u  is  governed  by  the  equation 

v2u  +  1  vv  .  u  =  0  ( 2 ) 

where  v  is  Poisson's  ratio.  The  gradient  and  Laplacian 
operators  in  three-dimensions  are  denoted  by  v  and  ,  re¬ 
spectively.  For  problems  exhibiting  symmetry  about  the  xy- 
plane,  which  contains  the  surface  of  discontinuity,  the  dis¬ 
placement  vector  u  may  be  expressed  in  terms  of  a  vector 
potential  $  with  components  $x,  <t>z  and  a  scalar  potential 

*  [5]: 

u  =  $  +  z  v  ip  (3) 

Hence,  it  is  not  difficult  to  verify  that  eq.  {2}  can  be 
satisfied  by  taking 


1 

3-4  v 


V ,  <f> 


and 


(4) 


V 2  <f>  =  0,  724»  =  0 
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The  displacement  vectors  for  problems  possessing  symmetry 
with  respect  to  the  yz-  and  zx-  planes  may  be  obtained  from 
eqs.  (3)  and  (4)  by  cyclic  permutation  of  the  variables 
x,y,z.  For  instance,  the  representation 


u 


+  XVi|)  ' 


3  ^  1 
3X 


1 

3-4  v 


V.  4>  ' 


(5) 


applies  to  problems  with  symmetry  about  the  yz-plane.  In 
eq.  (5),  <j> '  and  '  satisfy  the  Laplace  equation  in  three- 
dimensions. 


It  should  be  mentioned  that  eq .  (3)  or  eq.  (5)  is  a 
special  representation  of  the  more  general  solution  of 
Papkovitch  [6]: 


u  =  4 ( 1  - v )  B  -  V  (R  .  B  +  Bq) 


(6) 


where  R  is  the  position  vector.  Denoting  the  components  of 

B  by  B  ,  B  ,  B  ,  the  Papkovitch  functions  are  related  to  $ 

-  x  y  z 

and  in  eq .  (3)  as 


♦ 


x 


iio 

3  X  ’ 


(3-4v)Bz, 


t  = 


B 


z 


and  the  two  components  B  ,  B  are  taken  to  be  zero. 

x  y 

Once  the  displacements  are  known,  the  stress  tensor  o 
follows  directly  from  the  stress-displacement  relation 


=  u  [yjg—  (v.u)  I  +  vu  +  uv] 
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O 


(7) 


in  which  p  is  the  shear  modulus  of  the  material  and  I  is  the 
isotropic  tensor. 

Triaxial  Tension  Of  Elliptical  Disk 

Consider  an  infinite  solid  with  an  elliptical  disk  lying 
in  the  xy-plane.  The  z-axis  pierces  through  the  center  of 
the  disk  whose  surfaces  are  subjected  to  the  displacements 

Eux  =  -  [  o.|  -  v  { o2+03 )  ]x  » 

EUy  =  -  [ eg- v ( 03  +  0-1  )]y ,  Eu£  =  0  (8) 

for 

z  =  0  and  x2/a2  +  y 2  /  b  *  <_  1 

The  Young's  modulus  is  denoted  by  E.  Now,  the  negative  of 
the  displacements  in  eq.  (8)  correspond  precisely  to  those 
of  a  uniform  state  of  stress  in  a  solid  with  the  disk  absent, 
i  . e.  , 

°xx  =  °1*  °yy  =  °2 *  °zz  =  °3*  Txy  "  Tyz  Tzx  "  0 

Superposition  of  the  solutions  of  the  two  preceeding  problems 
will  leave  both  faces  of  the  disk  free  from  displacement  and 
will  yield  the  result  to  the  problem  of  a  thin  rigid  ellip¬ 
tical  disk  in  an  otherwise  uniform  state  of  stress.  Hence, 
it  suffices  to  solve  the  non-trivial  second  fundamental  prob¬ 
lem  owing  to  the  boundary  conditions  given  by  eq .  (8). 


-6- 


Let  f(x»y,z)  be  a  harmonic  function  such  that 


3  Z 


x  =  { 3 - 4 v )  =  ( 3-4 v )  |“  »  <f>z  -  0,  i|<  - 

From  eq.  (3),  the  displacements  become 


(10) 


3F  3F  _  32f 

ux  3 x  *  uy  =  3y  ’  uz  2  Jz7 


in  which  F  is  defined  as 
F  =  ( 3-4v )  f  +  z  || 


(ID 


Upon  substitution  of  eq.  (11)  into  (7)  gives  the  stress  com¬ 
ponents 


xx 

“ 


zz 

V 


h, 

2  u 


xz 


2u 


3  2  f 

3  2  F 

=  -  2v  ^ 

+  3  X2 

3  2  f 

3  2  F 

=  -  2(2-v) 

az2  + 

3  2  2  * 

3  2  f 

32  F 

=  -  2  ( 1  -  v ) 

3y  a  z 

3y3z 

3  2  f 

32F 

=  -  2 ( 1  - v) 

3  X  3  Z 

3  X  3  Z 

3  2  f  32F 


3  Z 


3y 


3  2  F 

3 x 3y  * 


(12) 


To  determine  the  only  unknown  function  f(x,y,z),  ellip¬ 
soidal  coordinates  C»  n,  c  will  be  employed.  The  rectangular 
coordinates  x,y»z  of  any  point  will  be  expressed  in  terms  of 
the  triply  orthogonal  system  5,  n>  t  in  the  form  [7] 
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a2 ( a2-b2  )x2  =  (a2+?)(a2+n)(a2+0 


b2(b2-a2)y2  =  ( b2  +  5) ( b2+n) ( b2+c ) 


a  2  b 2  z  2  =  5n5 


(13) 


where 


<»>5>0>n>  -  b 2  >_  5  -a 2 

In  the  plane  z  =  0,  the  inside  of  the  ellipse  x2/a2  +  y2/b2 
=  1  is  given  by  5  =  0,  and  the  outside  by  n  *  0. 

Making  use  of  eqs.  (11)  and  (13),  the  boundary  conditions, 
eq.  (8),  become 


(  3  -  4  v )  =  -  |=-  [o1-v(o2  +  o3)]x,  5-0 


(14) 


(3-4v)  "jw  ~  ~  £  to2~v^03+ol  ^  ® 


which  implies  that 


a  2  f 

=  constant,  5  =  0 

3  Z  2 


3  2f  +  3 2 f 

2  v  ^ 


The  solution  of  this  problem  can  be  obtained  from  the  known 
result  for  the  gravitational  potential  at  an  external  point 
of  a  uniform  elliptical  plate  [8],  i.e.. 


(15) 
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where 


Q(s)  =  s(a2  +  s) ( b  2  +  s  ) 

For  subsequent  use,  the  following  partial  derivatives  are 
computed : 

3  f  2A, 

Tx  =  J7^7  '  E(u)]x 

3 f  2A,  snu  cnu 

37  =  a^TT^  [E(u)  -  k’2u  -  k2  •  dn  u  ]  y  {16) 

The  variable  u  is  related  to  the  ellipsoidal  coordinate  £ 

by 

S  =  a2  ( sn~^  u-1) 

and 

u 

E  ( u  )  =  /  dn2 1  dt 
o 

The  quantities  snu,  cnu,  ---,  represent  the  Jacobian  elliptic 
functions  and  k,  k'  stand  for 

ak  =  (a2-b2)^^,  ak'  =  b 

A  glance  at  eqs.  (14)  and  (16)  shows  that  the  constants  A^  in 
eq.  (15)  cannot  be  evaluated  uniquely.  For  this  reason,  the 
additional  solution 

ux  =  -  A2x,  uy  =  A2y,  uz  =  0  (17) 
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will  be  introduced.  The  sum  of  eqs.  (14)  and  (17)  renders  a 
system  of  two  algebraic  equations  for  the  two  unknown  con¬ 
stants  and  A  2  which  yields 


ab2  ( 1  - v }  (  o-| +Oj> )  -  2v<jg 
A1  =  *  ETkT  *  4U  ( 1  +  v )  (  3- 4  v  ) 

(18) 

o  n  -  a  r>  3  -  4  v  a  2 

A2  =  -  jrifz  [<1  +  b^)E(k)  -  2K( k ) ]  A, 


where  K(k)  and  E { k )  are  the  complete  elliptical  integrals  of 
the  first  and  second  kind  associated  with  the  modulus  k, 
respectively. 


When  the  stress  state 


xx 


-  2yA2>  oyy  =  2yA2,  ozz  -  txy - 0 


is  added  onto  eqs.  (12),  the  contact  stresses  for  s  =  0  may 
be  calculated4.  The  normal  stresses 


(“xx1 


5  =  0 


V..0J 


o,-o,  3  (1 -v)  (o.+0„)  -  2va, 

*  [•]  (-V1)  -  7  t  nJ)(3-4v) — -] 


{ 1  -V  )  ( cr-j  +a2  )  -  2 vo ^ 


( “ z z  1  ^ _ q  ■  (1‘2v)  •  [ - (Hv)(3--40 


(19) 


4The  higher  order  derivatives  of  the  function  f(x,y,z)  can  be 
found  in  a  paper  by  Kassir  and  Sih  [9]. 


are  found  to  be  independent  of  the  geometry  of  the  elliptical 

disk.  For  n  =  0,  i.e.,  outside  of  the  ellipse  x2/a2  +  y2/b' 

=  i  a  a  and  o  become  singular  on  the  edge  of  the 
*  xx  *  yy ’  zz 

disk.  Further,  the  stress  exerted  by  the  surrounding  material 
on  the  disk  in  the  z-direction  vanishes  if  the  material  is  in¬ 
compressible.  The  shear  stresses  on  the  disk  are  given  by 


Ky> 


5  =  0 


=  0 


(1-v)  (oi+Ort)  “  2  vo0 

<’xz>f=0  -  2<1-'>)b  [  (Hv)(3-«v)t(k )-36.2/»2].(|-»2/.»-y2/b2) 

(  1  -V  )  (  Oi +Oo )  “  2  VOo 

<*yl>_  =  2d-v)b  c  (HV)(3-4v)ETiry^3^b23-(1-x2/a2-y2/b2) 


5  =  0 


(20) 


While  both  xxz,  xyz  are  zero  for  n  =  0,  they  are  unbounded  on 
the  boundary  of  the  disk  for  5  =  0  as  shown  in  eq.  (20). 

In  the  limiting  case  of  a  =  b,  E  =  K  =  tt/2,  the  constants 
A 1  and  A2  in  eq.  (18)  take  the  forms 

a3  ( 1  - v )  (o-j  +  o2)  -  2vog  ^  q2~v^°1^°3^ 

A1  =  "  *  (l+v)(3-4v)  ’  A2  2p  ( 1  +  v ) 

and  eqs .  (19)  reduce  to  the  results  for  a  penny-shaped  disk 
given  by  Collins  [1].  The  shear  stresses  in  eq.  (20)  may  be 
combined  to  yield 
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4(1 -v)  .  [ 


r/a 


o 


rz 


± 


( 1  - v )  { c-]  +<*£ )  -  2vo 
( 1  - v )  (  3-4 v ) it 


/I  -  ( r/a  )2  * 


0<r<a,  z  =  0 


where  a  =0  for  r  >  a,  z  =  0.  The  plus  and  minus  signs  refer 
rz 

to  the  upper  and  lower  faces  of  the  disk,  respectively. 

Returning  to  the  problem  of  finding  the  stress  distribution 
in  an  infinite  solid  containing  a  thin  rigid  disk  under  tri- 
axial  tension  at  infinity,  it  is  necessary  to  express  the  con¬ 
stants  A1  and  A2,  explicitly,  in  terms  of  the  applied  stresses 
at  infinity 

°xx  *  °r*  °yy  =  02°>*  =  °3~ 

which  are  related  to  o-j ,  o2»  o-j  in  eq .  (18)  as 

o-j  =  o<|-2yA2,  o2  ~  a2  +  2,jA2»  °3  =  03  (21) 

Inserting  eq.  (21)  into  eq .  (18),  it  can  be  easily  shown  that 
a,-.  02**  °3^  cannot  be  prescribed  independently.  This  re¬ 
striction  can  be  illustrated  by  considering  two  special  cases 
as  foil ows : 


Case  (i)  o  1  =  a  2  =  0 


Let  the  stresses  at  infinity  be 

°XX  =  =  -  2pA2’  °yy  =  °2°°  =  2yA2’  °ZZ  =  °3 
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^  A 


Solving  for  A-j  and  gives 


2pA1 


ab2  v 

( 1 +v }  ( 3-4v ) E ( k )  °3 


(22} 


2pA2 


n+VTk7 


[2-k2-2k' 2 


co 


°3 


Case  { i i }  a2  =  0 


Another  possible  solution  can  be  obtained  by  specify- 


i  ng 


xx 


=  a  ■ 


-  a 


1  -  2yA2 


yy 


=  c2*  =  2pA2, 


Z  2 


=  0 . 


It  follows  that 


a  3  k  2 [ VOg”  “  ( 1 “ V  )°1  ] 

2yAl  =  2(3-4v)L(l-v)K(k)  -  (l-va^/b^)E(Tn 


v  »  ®  n+v) 

2 v< A ^  =  o  2  =  2  ( o i  "*"a3  )  ”  2 


(23) 


[{a2/b2)E(k>  -  K(k)][vo3"-(l-v)o1“] 

*  - ( 1  - y  )  K ( k )  -  ( 1  - va2/ b2 ) E ( k } 

Eqs.  (22)  and  (23)  indicate  that  the  specification  of  the 
applied  stresses  is  severely  restricted  .  In  the  present  method 


5Such  a  restriction  was  also  mentioned  briefly  by  Eshelby  [4] 
in  his  survey  article  on  the  problem  of  the  ellipsoidal  inclu¬ 
sion.  3_ 


of  analysis  of  inclusion  problems,  it  appears  that  only  two 
of  the  three  principal  stresses  at  infinity  can  be  specified 
i ndependentl y . 

Elliptical  Disk  Displaced  Along  Its  Major  Axis 

Let  an  elliptical  disk  be  embedded  in  an  infinite  solid 
and  be  placed  in  the  xy-plane.  The  disk  is  displaced  along 
its  major  axis  by  the  amount  uQ,  a  constant.  The  necessary 
boundary  conditions  are 


u 


x 


,  t  =  0 


u 


z 


xz 


yz 


0  , 


n  =  0 


(24) 


The  symmetry  conditions  suggest  the  following  selection  of 
potential  functions: 


{ 3-4 v )  g  + 


ah 

ax ' 


*'  =  ttt 


ah 
ay  ’ 


ah 

7z*  * 


=  g 


(25) 


where  <j> '  ,  $ 1  ,  $ 1  are  the  rectangular  components  of  the  vector 
x  y  z 

<j>'  in  eq.  (5).  The  functions  g(x,y,z)  and  h(x,y,z)  satisfy 
the  Laplace  equations 

V2g(x,y,z)  =  0  ,  v2h(x,y,z )  =  0 

Putting  eq .  (25)  into  (5),  it  is  found  that 


4(l-v)g  +  if  ,  uy  =  ||  ,  uz 


14- 


3G 

az 


(26) 


From  eq.  (7),  the  components  of  stress  are  obtained: 


ag  a2G 

-2(2-v)^  + 


3g  3  2G 
lx  +  Ty7  ’ 


o  3g  32G 

TT  =  *2vax  +  Tz7  * 


3g  32G  a2G 

-2(l-v)j^  +  3x3y  »  2y  "  TyTz  * 


3g  3 2  G 

-  2  ( "I  »  v ) -yj  +  3X32 


(27) 


The  appropriate  harmonic  functions  for  this  problem  may 
be  chosen  as 


°»  d  s  2  B  i 

9(x,y’z)  *  B’  l 


(28) 


h(xty,z) 


ds 

(a2+s) /qTs7 


a  3k2 


[u  -  E(u)]x 


Note  that  h(x,y,z),  except  for  the  multiplying  constant, 
represents  the  derivative  of  the  gravitational  potential  at 
an  external  point  of  an  elliptical  disk  with  respect  to  x. 
For  the  purpose  of  evaluating  the  constants  and  B2,  the 
displacement  component  uz  is  computed: 


2x[nda2  +  S  )  (b2  +  5  }  ]  B2  n 

uz  - - abU-nJU^TJ  [  1  +  iTT^J 

The  condition  that  u  vanishes  everywhere  on  the  plane  z  =  0 


yields 


B  2  =  -a2B1 


(29) 


By  virtue  of  eqs.  (24),  (26)  and  (29)  for  «  -  0,  B,  is  found 

ak2  f30l 

=  "  2~  *  [  (  3-4\> )  k2  +  l  JK  (  k )  -  tO(T 


Knowing  B]  and  B2,  the  displacements  and  stresses  at  any 
point  of  the  solid  can  be  calculated.  On  the  plane  z  -  0, 
the  non- vani shi ng  displacements  are 


2  B  f 

<ux)  -  -  jp-  | [l+(3-4v)k2]u  -  E(u) 
a  (  kx)2  /  l  (  b 2 nT  . 

U-O  (az+iT  \l  az+^  y 

2B1xy  /  l  ’ 

(O  „  - - *  V  (a2H)  (bz  +  U 

J  n  =  0 


and  the  stresses  are 

8p  ( 1  -  v )  B  -j 


(  t  v  , ) 


xz 


C  =  0 


a  b 


( 1  -x2/a2-y2/ b2 ) 


(czz) 


1=0 


4p  ( 1  -  2  v )  B  i  x 

Ul  *  * 


b2  +  Z 

rta^nT 


(31  ) 


(32) 
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Both  t  and  o  are  singular  on  the  border  of  the  ellipse 

XX  Z  Z 

x2/a2+y2/b2  =  1,  while  t  =  0  everywhere  on  the  plane  z  =  0 . 

y  ^ 

When  a  =  b ,  K  =  E  =  x/2,  eq.  (30)  simplifies  to  the  form 
2au„ 

r  -  o 

B1  "  it (7 -8v) 

It  can  be  verified  that  for  r  >  a,  z  =  0,  s  ■*  r 2  -  a2,  and 
u  ->  sin'* 1  {-),  eqs .  (31)  and  (32)  are  in  agreement  with  eqs. 
(23)  and  (24)  in  [2],  respectively,  except  for 


( o 


zz 


) 

z  =  0 


8y(l-2v)  u 

- t— ) 

it  (  7  -  8v )  a 


cos  e 

(r/a)/(r/a)^ 


r  >  a  ( 33 ) 


where  uQ  corresponds  to  A  in  [2]. 


6Eq.  (33)  may  also  be  derived  directly  from  eq.  (20)  in  [2] 
if  the  order  of  integration  and  differentiation  is  properly 
observed  as  follows: 

1  3  +a  f(t)  dt  8yu 

(o  )  =  — ( 1  -2v ) —  [lim  /  -  -  •]»  f(t)  =  -  — - - 

zz  z  =  0  2  3x  z->-0  -a  /r*  +  ( z+i t T2  tt(7-8v) 

Carrying  out  the  integration  gives 


(°2Z)  = 
zz  z  =  0 


8u(l-2v)u  3 

71  (  7-8vl  3X  *-Sin 


(?)] 


8p ( 1 -2v )uQ 

=  7t(?-M 


(?)(r2-a2) 


1/2 


cos 


e 


Hence,  the  factor  (1-v)  in  eq.  (24)  of  [2]  should  be  replaced 
by  cos  e. 


-17- 


The  foregoing  method  of  solution  may  also  be  used  to 
solve  the  problem  of  an  elliptical  disk  displaced  in  an 
arbitrary  direction  by  a  constant  amount,  say  <5Q.  If  w  de¬ 
notes  the  angle  between  the  x-axis  and  the  direction  along 
which  the  disk  is  caused  to  move,  then  the  boundary  condi¬ 
tions,  eq.  (24),  may  be  generalized: 


ux  =  6o  cos  “*  uy  "  6o  s1n  <1)*  uz  =  °»  £  =  0 


u 


z 


n  =  0 


The  displacements  are  expressible  in  terms  of  four  harmonic 
functions  as 


u 


x 


-4(1 -v)g1 


+ 


I!° 

ax  ’ 


u 


y 


■  4  ( 1  -  v )  g  2  + 


!!< 

ay 


!!o 

9  Z 


in  which 


Go  =  61  +  G2»  G1  =  X91  +  hl*  and  G2  =  yg2  +  h2 
To  satisfy  the  Laplace  equations  in  three  dimensions,  g.(x,y,z) 

J 

and  h.(x,y,z)  are  taken  in  the  forms 

J 


<  r  7  ds 

g.(x,y,z)  =  C.  /  - 

J  J  C  /QlsT 


,  j  =  1,  2 


h^x.y.z)  =  D-j  x  / 


ds 


£  /a^+s  7qT?T 


h?(x,y,z)  =  D?y  / 

e 


_ ds 

( b2+s ) /Q { s ) 
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Since  the  displacement  uz  vanishes  for  z  =  0 ,  the  constants 


D1  =  -  a 2  C  -j  ,  D2  =  -b2C2 

The  remaining  unknowns,  say  Cj  (j  =  1,  2),  can  be  evaluated 
from  the  boundary  conditions  yet  to  be  satisfied  and  the  solu¬ 
tion  of  the  problem  is  essentially  complete. 

Displacement  Of  Rigid  Sheet  With  Elliptical  Hole 

Suppose  that  two  semi-infinite  solids  are  bonded  perfectly 

to  a  thin  rigid  sheet  with  an  elliptical  opening  through  which 

the  solids  are  connected.  The  sheet  is  allowed  to  move  in  the 

plane  z  =  0  by  a  constant  amount  parallel  to  the  x-axis.  The 

equivalent  condition  is  to  specify  a  constant  shear  stress 

t  =  x  for  e  =  0.  For  this  problem,  the  following  conditions 
ZX  0 

must  be  satisfied: 

u  =  uy  =  0,  n=0;  u _  =  0,  z  =  0 
x  y  *■ 


(34) 


T 


The  problem  may  be  formulated  in  terms  of  a  single  function 
p(x,y,z)  which  is  related  to  <b  and  in  eqs.  (3)  and  (4)  as 


*x  =  *y 


where 


v2p(x,y,z)  =  0 


The  representation  of  the  components  of  displacement  as  given 
by  Trefftz  [5]  is 


UX  =  +  z  uy  *  2  ~  ^ 


,*,y  ■  2  HIT  1351 


The  stresses  corresponding  to  eq.  (35)  are  given  by 


o  3  3  P  82P  o  3  3p  32P 

=  Tx  t-(3'2v^TI  +  2  Tx7^  ’  =  ?x  ^'2v3z  +  2  sy7-*  ’ 


o  3  3  p  3 2  P  T  3  3p  P 

#  ■  37  Cd-Zv)^  +  z  -P  ■  57  +  22 


3  x  3y  ■ 


t  ,  32  3  p  t7¥  3 2  P 

y  —  =  -  [  p+2z  - — ]  ,  -  -  -  ( 3 - 4  v ) 

u  3x3y  LM  3 z J  v 


3  2  3  p 

3^  +  377  Cp+2z  37] 


(36) 


On  the  plane  z  =  0,  eq.  (34)  requires  that 


12.  = 
3  Z 


0, 


n  =  0 


( 3-4 v ) 


e  =  0 


(37) 


The  first  condition  in  eqs.  (37)  is  satisfied  automatically  by 
taki ng 


p(x,y,z)  ■  y~  \  tr? 

5 


+  s 


y2 

PT? 


i] 


ds 

swm 
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while  the  second  condition  yields 


a3k2k ' 2  t  Q 

2uC  =  k,2K(k)+[(3-4v)k^-k'^]E(k) 

Once  p(x,y,z)  is  determined,  the  displacements  and  stresses 
throughout  the  solid  can  be  computed  from  eqs.  (35)  and  (36). 


For  z  =  0,  both  uy  and  uz  vanish  and 


(uj 


5  =  0 


2C  <3~-4-v-*  (l-x2/a2-y2/b2)1^2,  (uv)  =  0 


‘n=0 


The  stresses  on  the  plane  z  =  0  are 


(o„)  =  -  ( 1  -x2/a2-y2/b2  )~  1/2 


zz 


5  =  0 


(ty2) 

yz  n  =  0 


2uC  xy 


(t,v)  =  2yC 


(5-c)«/QTTT 

3-4 v  r  ab2 


zx 


n  =  0 


{■ 


[ 


ab2  /QTZT 


E(u)  + 


snu  cnu 

dn  u 


+  u - E ( u )  _  x; 


a J  k  2  ‘  U-c)(a 


I  b  2n^ 

7+Tl  V  5(a2+5) 


and 


(38) 


(  °zz  )  =  ^TVZ^  "  °*  ^zx^ 

zz  n=0  yz  c=0  zx 


=  T 


Using  L'  Hospital's  rule,  the  constant  C  for  a  circular 

t 

hole,  a  =  b,  may  be  recovered: 

2  3  3  T 

r  =  _ 0 

u  hu(7-8v) 
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Aside  from  a  couple  of  misprints,  (u  )  ,  (t  )  ,  and 

x  5  =  0  y  n  =  0 

(T  )  check  with  those  given  by  eqs.  (41)  and  (42)  in  [2] 
zx  n=0 

if  t  is  identified  with  o  .  The  expression  for 
o  o 


<°zz> 


_  _  8^  1  - 2\.?.  — -  T ^  cos  9 

z=0  u(7-8v)  /l - ( r/a ) "  0 


fails  to  agree  with  that  of  [2]  for  the  same  reason  as  men 
tioned  earlier  in  footnote  (6). 


Axial  Displacement  Of  Elliptical  Disk 

If  a  thin  rigid  disk  of  elliptical  shape  is  given  a  con 
stant  displacement  wQ  normal  to  its  plane,  then 


ux  =  uy  =  °,  z  =  0;  uz  =  w0,  5  =  Q 


(39) 


which  suggests  that 


’x  ”  *y  =  °’ 


- ( 3-4v ) q ,  4*  =  q 
Inserting  eq.  (40)  into  (3),  the  result  is 


(40) 


ux  ~  z  lx’  uy  =  z  ay*  uz 


( 3-4v)q  +  z  ff" 


(41) 


From  eq,  (7),  it  is  further  found  that 


a2q 


3q  a2q  °vv  3q 

V"  ="2“Tz  *  z  TP-  2^  *  •2',sl  +  z  W 


°zz  3  q  ^  q  T x~y  -  7 

IT  =  -2(1-v)3i +  z  —  ■  ^  -  2 


JT7'  2y 


a2q 

3X3y 
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3q 


1  1/  7  ^ 

dr  - 


32q  t 


3q 


a2q 


ay  +  2  ayaz’  2p  ‘^'2v^3x  +  z  3xaz 


The  only  unknown  function  q{x,y,z)  satisfying 


v2q(x,y,z)  =  0 


can  be  taken  in  the  form 


q ( x *y ,z )  =  D  J 


ds  2D 


(43) 


£  /Q ( s  j  a 

Eqs.  (39),  (41)  and  (43)  may  be  combined  to  give 
D  = 


aw 


2(3-4v)K(Y7 

Calculating  for  the  derivatives  of  q(x,y,z)  with  respect  to 
x,y  ,z  ,  i .e. , 


34  aw0x  / 

37  =  (3-4v)U-n)U-c)K{k)  *  ^ 


Ub2  +  S) 
a^+l 


3  4 


aV 


37  =  (3-4v)(e-n)(e-c)K(kT 


t(a2+c) 

b^n 


,  ^  1  /  2 

7z  =  ( 3-4v ) b  (5-n)(5-t)K(k)  *  /U<t  +  U^  +  ?T 
and  so  on  the  non-trivial  displacements  and  stresses  for 

z  =  0  are 

w. 


(u*»t.0  ’  (“2>n=0  =  ^ 
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and 


‘ °zz 1 z  =  0  * 


=  + 


4U  C 1  - v ) wQ 
(3-4v)b  K(kT 


(l-x2/a2-y2/b2) 


-  1/2 


5  =  0 


“  "I 

(Txz)  + 
xz  z  =  0 

2u(1-2v)wq 

/( a2  +  ; ) ( bz  +  i  ) 

[ - 

H 

Nf 

Nl 

11 

O 

+ 

_ 

(3-4v)51/2(^-i;)k  K(k) 

/(  a2  +  c  )  L-  (  bz  +  c  )  J_ 

(44) 


0 


in  which  -(b2n)  is  a  positive  definite  quantity.  The  nota¬ 
tions  z  =  0+  and  z=0~  refer  to  the  upper  and  lower  faces  of  the 
disk,  respectively. 


The  force  exerted  by  the  elastic  solid  to  oppose  the 
displacement  of  the  elliptical  disk  may  be  found  from  the 
i nteg  ra 1 


F 

z 


•  U  U. 


z  =  0 


+ 


(o  )  ]  dxdy 

zz  z=0' 


The  region  £  is  bounded  by  the  ellipse  x2/a2+y2/b2 
stituting  eq.  (44)  into  (45),  Fz  is  obtained: 


(45) 

1.  Sub- 


F 


z 


8p { 1  - v  )  w  _w? 

( 3-4 v ) b  K(kT  V  n-Xz/az-yz/bM  1/2dxdy 

1 6 -n p  ( 1  - v )  awQ 
(  3  -  4  v )  K  (  k  ) 


(46) 


In  the  limit  as  a  ■>  b,  eq.  (46)  reduces  to  Collin's  solution 
[1]  for  a  circular  disk. 
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T  h  ree- D imensional  Stresses  Near  Inclusion  Border 


For  the  purpose  of  establishing  possible  failure  criteria, 
the  stresses  near  the  border  of  a  plate-like  inclusion  will 
be  investigated.  It  is  convenient  to  introduce  a  rectangular 
cartesian  coordinate  system  n,t,z  such  that  the  origin  of  this 
system  traverses  the  periphery  of  the  inclusion.  The  zn-, 
nt-,  and  tz-  planes  are  known,  respect i ve 1 y ,  as  the  normal, 
rectifying  and  osculating  planes  to  the  curve  which  will  be 
taken  in  the  form  of  an  ellipse. 

In  the  immediate  vicinity  of  the  inclusion  border,  the 
ellipsoidal  coordinates  £ ,  n,  ?  can  be  expressed  in  terms  of 
the  polar  coordinates  r,  0  defined  in  the  nz-plane,  where  r 
is  the  radial  distance  measured  from  the  edge  of  the  inclusion 
and  0  is  the  angle  between  r  and  the  n-axis.  The  required  re¬ 
lationships  of  n»  ?  to  r,  e  are7 


e  = 


2abr 


(a2sin2<(>  +  b2cos2<f>)"'^ 


cos 


2  2 


2abr 


n  =  - 


( a  2  s i n2  $  +  b2cos  2  $ ) ^  72 


s  i  n 


2  2. 


C  =  -  (  a  2  s  i  n  2  4>  +  b  2  c  o  s 2  $  ) 


(47) 


7A  detailed  derivation  of  eq.  (47)  is  given  in  [9]. 
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In  eq.  (47),  r  is  assumed  to  be  small  in  comparison  with  a 
(or  b)  and  <j>  is  the  angle  appearing  in  the  parametric  equations 
of  the  ellipse,  i.e., 

x  =  a  cos  y  =  b  sin  ^ 

Since  the  derivation  of  the  local  stresses  is  similar  to  those 
given  by  Kassir  and  Sih  [9]  for  the  three-dimensional  crack 
problem,  the  detail  calculations  will  be  omitted  here.  By 
means  of  eq.  (47)  and  the  appropriate  equations  for  finding 
the  stresses,  the  following  results  are  obtained: 


k,  e  e  3e 

=  +  —  cos  -  ( 3- 2 v  -  sin  -  sin  — ) 

/2r  2  2  2 


k?  e  e  3e 

+  —  sin  —  (2v  +  cos  —  cos  — )  +  0(1) 

/27  2  2  2 


k,  e  e  3e 

=  -  —  cos  -  ( 1  - 2 v  -  sin  -  sin  — ) 

/2r  2  22 


k  0  0  3  0 

+  —  sin  -  ( 2- 2 v  -  cos  -  cos  — )  +  0(1) 

/2?  2  2  2 


=  +  —  .  2v  cos  -  +  —  2v  sin  -  +  0(1) 

/2?  2/2?  2 


T 


nt 


—  cos  -  +  0(1) 
/2r  2 
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T 


n  2 


k,  0  0  30 

+  —  sin  -  ( 2-  2v  +  cos  -  cos  — ) 

/If  2  22 


k9  e  0  3e  ,  v 

+  —  cos  —  ( 1  - 2 v  +  sin  -  sin  — )  +  0 { 1 ) 
/2r  2  2  2 


k  0 

=  _2  sin  -  +  0(1)  (48) 

tz  /27  2 

Although  these  stresses  were  derived  from  the  solution  of  an 
el  1 iptical ly-shaped  inclusion,  they  are  in  general  valid  for 
a  plane  inclusion  of  arbitrary  shape.  Moreover,  the  inclusion- 
border  stress  fields  for  the  four  preceding  boundary-value 
problems  are  included  in  eq.  (48)  as  special  cases. 

Now,  it  is  significant  to  observe  that  eq.  (48)  is  com¬ 
posed  of  the  linear  sum  of  three  distinct  stress  fields  each 
of  which  can  be  associated  with  a  different  mode  of  deformation. 
Referring  to  Figs.  1(a)  through  1(c),  the  intensity  of  the  local 
stresses  at  the  point  P  caused  by  the  movements  of  the  inclusion 
in  the  n- ,  z - ,  and  t-  directions  are  governed,  respectively,  by 
the  three  parameters  k-|  ,  k^  and  k^-  These  three  modes  of  dis¬ 
placements  are  necessary  and  sufficient  to  describe  all  the 
possible  displacements  of  the  inclusion.  It  will  be  shown  sub¬ 
sequently  that  the  parameters  kj  ( j  =  1 ,  2,  3)  depend  only  upon 
the  prescribed  stresses  or  displacements  and  the  inclusion  ge¬ 
ometry.  The  singular  behavior  of  the  i nc 1  us i on- border  stresses 
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is  the  same  as  that  for  a  sharp  crack.  In  other  words,  the 
1 / /r  type  of  stress  singularity  is  preserved.  However,  un¬ 
like  the  crack  problem,  the  angular  distribution  of  the 
stresses  is  a  function  of  the  Poisson's  ratio  of  the  elastic 
sol  id. 

A  close  examination  of  the  stress  expressions  in  eq.  (48) 

reveals  that  a  ,  o  ,  and  t  correspond  precisely  to  those 
nn  z  z  n  z 

obtained  by  Sih  [ 1 0 ] 8  for  a  rigid  line  inclusion  under  the 

conditions  of  plane  strain*  In  fact,  the  stress  component 

o  is  equal  to  v  ( a  +  a  ) ,  a  condition  which  is  well  known 
1 1  n  n  z  z 

in  the  analysis  of  plane  strain  problems.  The  shear  stresses 

i  and  -t.  can  be  identified  with  the  two-dimensional  problem 
n  t  tz 

of  a  line  inclusion  subjected  to  longitudinal  or  out-of-plane 
shear  loads.  Hence,  the  stress  state  around  a  plane  inclusion 
in  three-dimensions  is  locally  one  of  plane  strain  combined 
with  longitudinal  shear. 


8The  stresses  orr,  c  ,  and  ir0  given  by  eq.  (48)  in  [10] 

should  be  transformed  into  rectangular  components  oxx,  a  , 

t  in  accordance  with 
xy 


o  +c 

xx  yy  rr  ee 


°yy  ■  °xx  +  21’xy  =  e'2’  ° +  2iT-) 


ee  rr 


re 


For  K  =  3-4v,  the  functional  forms  of  o  ,  o  ,  t  correspond 

to  a  ,o  ,  -r  in  this  paper,  respect  i  V  e  1  yy.y  y 

nn  z  z  nz 
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In  general,  the  three  parameters  k ^  (j  =  1,  2,  3)  will 
occur  simultaneously  over  the  inclusion  border.  They  may  be 
i nterpreta ted  as  a  measure  of  the  elevation  of  stresses  due 
to  the  presence  of  thin  rigid  inclusions  embedded  in  elastic 
solids.  From  eq.  (48),  the  formulas 


k 


1 

1  = 


1  i  m  /2 "r  (  a 
r-0  z 


e=o 


k 


2 


1 


1  im  /2 7  ( t 
r+0 


e  =  0 


(49) 


k,  =  1  im  /2r  (t +  ) 

J  r+0  1  9=0 

are  obtained.  Eq.  (49)  may  be  applied  to  evaluate  kj  for  the 
bo u ndary - va 1 u e  problems  solved  earlier.  Following  the  work 
of  Kassir  and  Sih  [9],  it  is  found  that 


(1)  Triaxial  Tension. 


=  + 


( 1  -v )  ( o-j  +^2 )  -  2  voj  b 
(1+v) (3-4v)Elkj“ 


1/2 


(a2sin2<j>  +  b2cos2<)> ) 1  , 


k2  k3  0 


(50) 


(2)  Parallel  Displacement. 

2pak2u  b 

kl  ■  '  L(3-4v)k2-MjKtk)-b(k)  (I>  U2si"J* 

+  b2cos2,)”^a  cos  k.,  =  0 
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4u ( 1 -v)ak2u 


‘3  =  [(3-4v>)k2  +  lJK(k)-E  (~kT 


1/2 

(tt)  ( a2 s i  n2 4 


+  b2cos2$)“  3^4  sin 


(51  } 


(3)  Rigid  Sheet 


2  b  k  2 


1/2 


kl  "  + 


[(5-4v)k2-k''2]E(k)  +  k'7l<Tk7  (T’  (a2sin'!* 


+  b2cos2$)"  1/121  cos  * ,  k9  =  0 


(3-4v )ak2T 


;3  "  L  (3-4v)k^-k' ^ JE( k)  +  k ,  ZK( M  ^a 


1/2 

(-)  (a2sin 


+  b2cos24> )  ^  4  sin  $ 


(52) 


(4)  Axial  Displacement 


2vtw 


1/2 


k-|  -  0 ,  V>2 


=  •  (  3  -  4  v  )  K  (TcT  (b>  (a2sln2$ 


+  b2  cos2  <j> )  ~  ^  4 ,  k,  =  0 


(53) 


It  is  interesting  to  note  that  k,  are  not  constants  but  func- 

J 

tions  of  position.  E q .  (50)  is  associated  with  the  local  dis 
placement  shown  in  Fig.  1(a)  while  eq.  (53)  with  Fig.  1(b). 
The  displacement  modes  pertaining  to  the  results  in  eqs.  (51) 
and  (52)  are  more  complicated.  For  0  <  4  <  inclusion 
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border  experiences  a  combination  of  the  movements  illustrated 
in  Figs.  1(a)  and  1(c).  The  parameters  k-j  and  k3  attain  their 
maximum  values  at  $  =  0  and  respectively. 

For  problems  involving  all  three  parameters  k^  (j  =  1,  2, 
3),  it  is  possible  to  postulate  a  criterion  of  failure  for 
rigid  inlcusions  in  the  form 

fcr  =  f (k1 ,  k2,  k3) 

which  states  that  failure  of  the  material  surrounding  the  in¬ 
clusion  occurs  when  the  combination  of  k^  ,  k2,  and  k3  attains 
some  critical  value. 
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